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ABSTRACT 


We show that the electron streams that give rise to 
Type III solar radio bursts are stable and will not be decele- 
rated while propagating out of the solar corona. The stabili- 
zation mechanism depends on the parametric oscillating two-stream 
instability. Radiation is produced near the fundamental and second 
harmonic of the local electron plasma frequency. Estimates of the 
emission at the second harmonic indicate that the wave spectra 
created by the oscillating two-stream instability can account for 
the observed intensities of Type III bursts. 



I. Introduction 


Sporadic Type III solar radio bursts are caused by particle 
streams accelerated in the vicinity of a chromospheric flare. These 
energetic particles propagate through the solar corona at a constant 
velocity between 0.3 - 0.6 c, and excite radiation at the fundamental 
and second harmonic of the local electron plasma frequency, . 

Recent observations at 1 AU using data from the radiometer and particle 
experiments on the IMP-6 satellite (Alvarez, Haddock, and Lin 1972; 

Lin 1973) demonstrate that the particle streams consist primarily of 
electrons with energies in the range 1-100 keV. 

The interaction of such streams with the background coronal 
plasma produces unstable electrostatic waves with frequency near the 
electron plasma frequency and phase velocity near the beam velocity 
(Sturrock 1964; Kaplan and Tsytovich 1968; Zheleznyakov and Zaitsev 
1970b; Smith and Fung 1971). According to current views, these electro- 
static waves are subsequently transformed into electromagnetic radiation 
at U) e by scattering from ion density fluctuations, and at 210 ^ by 
scattering from each other. A review of the observations and theory can 
be found in Smith (1970) . 

While the above qualitative picture is generally accepted, 
quantitative estimates of the radiation spectrum and intensity require 
a more precise nonlinear theory of the interaction of the beam with the 
plasma. To estimate the electromagnetic radiation and the spectrum emitted 
in Type III bursts, a knowledge of the plasma wave energy spectrum before 
transformation is required. This in turn depends on whether the beam is 
decelerated to form a plateau in velocity space, or whether instead non- 
linear effects cut off the plateau formation at an early stage and sta- 


bilize the beam. 
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Kaplan and Tsytovich (1968) calculated the effect of induced 
nonlinear scattering of the plasma waves by thermal fluctuations, and 
concluded that this scattering was fast enough to stabilize the beam 
before the formation of a plateau. However, as Zheleznyakov and Zaitsev 
(1970a) correctly pointed out, the calculation of the scattering time 
by Kaplan and Tsytovich was inaccurate. Zheleznyakov and Zaitsev cal- 
culated the correct time for the development of this nonlinear process 
and concluded that it was longer than the characteristic time for quasi- 
linear deformation. Therefore, in contrast to the conclusion reached 
by Kaplan and Tsytovich, stabilization does not occur. In subsequent 
publications, Zheleznyakov and Zaitsev (1970b) and Zaitsev, Mityakov 
and Rapoport (1972) calculated the power emitted in Type III bursts, 
based on the hypothesis that nonlinear effects do not stabilize the beam. 

Smith and Fung (1971) also considered the question of whether 
beam-induced plasma waves could be scattered by polarization clouds of 
ions to stabilize the beam. They concluded that this mechanism could 
not stabilize electron streams in the solar corona. 

One of the most striking features of Type III bursts, however, 
is that the source region propagates at almost constant velocity from the 
inner corona to nearly 1 AU (Fainberg, Evans, and Stone 1972). In this 
paper we demonstrate that a nonlinear mechanism does exist that will pre- 
vent the deceleration of the electron stream. The mechanism that we 
consider is the "oscillating two-stream instability" (Nishikawa 1968a, b; 
Sanmartin 1970). This instability is driven by the fields generated 
through the beam-plasma instability, and transfers the energy in these 
fields to regions of phase space where the fields do not interact reso- 
nantly with the electron stream. The result is to prevent deceleration 
of the stream by its self-generated fields. 
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In the next section, we discuss the theory of the oscillating 
two-stream instability. In Section III, we calculate the conditions 
under which this nonlinear process can prevent plateau formation. In the 
last section, we use these results to estimate the power emitted in Type 
III bursts. However, we defer to a later paper a detailed consideration 
of the frequency spectrum and its temporal evolution. 


II. The Oscillating Two-stream Instability 

This instability has been discussed recently in the literature 
as an efficient mechanism for coupling energy from high-phase-velocity 
oscillations into shorter-wavelength oscillations which readily transfer 
energy to the ambient coronal plasma. A detailed mathematical treatment 
based on the fluid equations can be found in Nishikawa (1968a, b); an ana- 
lysis based on the Vlasov-Maxwell equations has been given by Sanmartin 
( 1970 ) . 

To our knowledge this instability has not yet been discussed 
with respect to astrophysical plasmas, so we present here a brief but 
self-contained derivation of the dispersion relation based on simple 
physical arguments and refer the mathematically interested reader to the 
original references. 

Consider a plasma in the presence of a long-wavelength high- 
frequency field 


E = E cos (w t + k x) 
o o o 

such that the phase velocity of the wave w 0 ^ 0 >> » the thermal 

velocity of the plasma. 

In the absence of this oscillating field, the normal electro- 


static mode of- the plasma is the electron plasma oscillation with 
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2.2% 

wavenumber k and angular frequency to ^ = (0^(1 + 3k A^) , and 

the ion acoustic wave with the same wavenumber and with frequency 
2 2k 

to. = kC / (1 + k to) 2 . (In the unperturbed plasma with equal ion and 
A s D 

electron temperatures , this mode is heavily damped and cannot properly 

be considered a normal mode.) Here the ion sound speed C g = [ (y^KT^ + 

k 2 % 

Yi KT i)/M] 2 ; the Debye length A^ = (KT e /47m e e ) ; e is the unit 

(positive) electronic charge; K is Boltzmann’s constant; n e > T » 

and T. are respectively the electron density, electron temperature and 

ion temperature of the coronal plasma; y and y^ are the adiabatic 

2 2 

indices of the electrons and ions; to = 4ime /m ; and m , M are 

e 

the electron and ion masses, respectively. These normal modes obey the 
equations 


3t _2 n ek + V 1 3t n ek + W ek n ek ° ’ 


( 1 ) 


9 9 2 

— 9 n .M + -r— n + to n = 0 

ik 2 St lk A xk 


where n ^ , n^ are the spatially-Fourier-analyzed electron and ion 
density fluctuations with wave number k , and are the total 

damping rates (collisional and collisionless) of these fluctuations 
respectively. The presence of an external force due to a long wave- 
length field oscillating at a frequency (0 q 'v to couples these two 
modes. For simplicity in calculating the coupling forces we assume 
k Q = 0 . (The extension to finite k Q is straightforward and does not 
influence our conclusions for the cases of interest here.) 

The coupling force can be derived by simple arguments. Let 
n^(x) , n g (x) be the ion and electron density fluctuations. The driving 


9 
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field tends to displace the electrons relative to the ions. This driver 

now produces high frequency charge fluctuations with a frequency 03^ . 

Then, because n - n. , we have 
e 1 


3n. 

6n = t ~ 
e 3x 


-3n. eE 


Ax = 


3x 


m03 


— COS 03 t 
2 o 


= ikn . , 
ik 


eE cos 0) t 
o o_ 

2 

m 03 


where we have used a Fourier space transform to replace 3/3x by ik . 
The electron equation now becomes 


3t 


2 n ek + V 1 3t n ek + W ek n ek 


ikn.. eE cos 03 t 
ik o o 


m 


( 2 ) 


The coupling force on the ions is due to the gradient of the electric 
field pressure. 


_1 3 _ 

8tt 3x 


(E cos 03 t+ E(x) 
o o 


y 


where E(x) is the self-consistent electric field. Using Poisson's 
equation, we have 


2 

~ (E cos 03 t + E(x) ) 
a o O 

3x 8ir 


-eE 6n cos 03 t 
o e o 


2 

The term E (x) is neglected. The equation for the ion acoustic fluc- 
tuations becomes 


3t 


2 n ik + V 2 


a , 2 1 

7T“ n., + 03. n.. = — 

3t ik A ik M 


3x 


[eE 6n (x) cos 03 t] 

n o n J 


ikn , eE cos 03 t 
ek o o 

M 


( 3 ) 
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The coupled equations (2) and (3) describe the instability. 
To derive the dispersion relation for the system, we Fourier transform 
equations (2) and (3) in time according to 


” e k (t> 


\ 


n ik (t) / 


doj e 


-iu)t 


n ek (“) 


n± k < U ) 


and obtain the system 


(a)2 ‘ w ek + iv l w) n ek (a)) + i <k [n ik ( ^o ) + n ik (a>Hi) o )] = ° 


(a) i 2 -w^ + iv 2 w) n ik (w)’ - ifi 2 k [n ek (a>-w o ) + n ek (uHu o )] = 0 


(4) 


where = keE Q /2 (i/ m , i/ M ) . We see that couples 

with n ek (w + u> q ) , and n ek (u> + C0 q ) couples with n ik (w) and 

n ik ( “ - 2<i) o ) * 

We then obtain the (infinite) system of equations 

r 

(w 2 -a) 2 + iv 2 w)n ik (w) - i«i k n ek (w-W o ) - ifl^ek^o 5 = 0 


ifi ek n ik (w) + 


4 + iv 1 < w - W o )]n ek (a,_a, o ) + ifi ek n ik (a, ~ 2w o ) 


= 0 


i 4 n ik (w) + ifi ek n ik (aH ‘ 2w o ) 


+ [(aH-w o ) 2 - w 2 k + iv 1 (a>fto o )]n ek ((vHo o ) = 0 


V. 


We shall be interested in Re w ^ ^ «tO Q and thus neglect n^aH^) 


( 5 ) 
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as being off-resonant. This closes the infinite set (5), and the disper- 
sion relation for the resonant normal modes of the coupled system is given 
by the secular equation 


ft) 2 - ft) 2 + iv ft) -iftj " ±Q ±k 

A i ik 

^ek (0>-«> o ) 2 -4«V“-“o> 0 


- 0 ( 6 ) 



0 


2 2 

(oj+ft) ) - to , +iv.. (ftrtto ) 

o ek 1 o 


Assuming that ft)^ = u)^ , and that the damping rates v ,v ^ are small, 
we may approximate 


(w 


+ 0 ) ) 2 - 
— o 


“ek + iu l<“ 


+ ft) ) * 
— o 


2 “ek ( “ 


+ 6 ) + 


ift) v. 

o 1 


where 6 = ft) Q - ft) ^ is the "frequency mismatch" between the driver and the 
electron normal mode. 

The secular equation (6) becomes simply 

k 2 e 2 E 2 

(CO 2 + iow 2 - CO 2 ) [(10 + iv ) 2 - 6 2 ] + -55^ jp_ - o . (7) 

ek 

We now write ft) = ft) + i V , where ft) and y are real. In the 
present analysis, y > 0 implies instability, which is reflected in the 
growth of ion acoustic waves at frequency ft) and electron plasma oscilla- 
tions at ft> r + ft) . From (7) we see that we may have y > 0 under the 
following conditions : 

(i) 6 < 0 , ft) f = 0 (i.e., y » ; 

(ii) 6 > 0 , ft) + 0 

r 
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Case (ii) is the resonant three-wave interaction of the high-frequency 
driver with two lower-frequency normal modes, known as the "decay in- 
stability" (Tsytovich 1970, Davidson 1972). We shall be concerned with 
case (i) , because in the case of Type III bursts the driver wave is simply 

the beam-excited plasma wave with cd - to . , and therefore 6 = id - cd , 
r o e o ek 

= 0 ) e - U) e k < 0 . From Eq. (7), with 6 < 0 , we outline below the basic 
features of the instability. 

There is a threshold amplitude for the instability to occur, 
which is found by setting cd^_ = y = 0 in Eq. (7) and minimizing with res- 
pect to 6 . The result is 


V. 


8tt n KT 
e e 


= 4 


CD 


( 8 ) 


ek 


2 2 

where we have used the fact that k L « 1 and have assumed that T - T. 

D ex 

Notice that represents the total electron damping and, for 


the case of a Maxwellian plasma, can be written as 



(9) 


where is the electron-ion collision frequency. 

The growth rate for fields well above threshold is given by 


lk 


°ek ( M ) 


1/3 


<^ 2 d> 


1/3 



( 10 ) 


where we have maximized y^ with respect to 6 and neglected 
and (D^ compared with y^ . 

Equation (10) describes the growth rate of the oscillating 
two-stream instability at very early times. We shall also be interested 
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in the radiation produced after the plasma has reached a quasi-steady 
state and is near marginal stability. The growth rate then becomes 


2 2 

k e 2 2 

Y lk = 4mM (E o “ E T ) 


( 2 “ ek ) C 2 _ v v 12 ) 

A 12 


(id 


The electron plasma waves produced as a result of the oscillating 

two-stream instability have much lower phase velocities than the driver , 

wave. This can be seen from the energy and momentum conservation conditions 

2 2 

for the waves, plus the usual restriction that k ^ « 1 . Denoting the 

driver, electron and ion waves by (oj q , k Q ) , ( w ek »k ), and ( w » 

the conservation relations are GJ - to + oj and k - k..+ k . Noting 

ore o i e 

that to = to, - (0 , k = V. /V , to. -to , k =* aA _ , 

o ek e o beD ek e’e D ’ 
o e 

to « Y 1t « to , k.=k - k - - k , with a - .1 - .2 and V, - 20 V , 
r'lk e xoe e b e 

e 

we find that the phase velocity of the growing wave is of order ~ 5-10 , 

while the phase velocity of the driver is w e /^ 0 - 20 V g . Notice that the 
excited ion waves appear only as solutions of the nonlinear dielectric 
constant and disappear in the absence of the driving field. 

Because of. the above interaction, the energy of the stream-excited 
plasma waves is transferred to a region of velocity space where these waves 
do not interact resonantly with the stream particles. We stress, however, 
that the energy transfer is not direct; the nonresonant waves are induced 
through the driven bulk motion of the electrons relative to the ions, and 
are not in random phase with the driver. Thus, the energy and momentum 
relations above need only be satisfied as approximations. In the non- 
resonant region of velocity space, the induced plasma waves can be Landau- 
damped by the ambient plasma, creating energetic tails in the ambient 
electron distribution. Some of the consequences of this tail formation 
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with respect to the Type III radiation will be examined in Section III 
below. 

Before closing this section, we point out that the above theory, 
including such consequences as tail formation, has been confirmed by a 
series of computer simulations following both electrons and ions in their 
self-consistent orbits (Kruer and Dawson 1970, 1972; Kainer, Dawson, and 
Coffey 1972) . 


III. Nonlinear Effects on the Beam-Plasma Interaction. 

We examine next the collisionless relaxation of an electron beam 
propagating in the solar corona under the following assumptions: 

(i) the interaction is one dimensional, i.e. all the plasma 
waves in the system propagate in the direction of the 
beam. 

(ii) the excitation of the plasma waves is considered under con- 
ditions that the beam particles can resonate with waves 

in a finite range of wave number k . This implies the 
usual restriction 


AV, 


(5 




1/3 


/ 


( 12 ) 


for the validity of a quasilinear description, where AV. 

b 

is the velocity spread in the beam and n^ is the density 
of the beam. 

(iii) The system is homogeneous. 

For a discussion of these assumptions, see Smith (1970). 

Under these assumptions the evolution of the beam will be described 
by the well-known set of quasilinear equations 
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3F b (v,t) 

3t“” 


■ c) b V v - c) ] 


(13) 


3& ok< t > 


9t 2 Y 2k^ok (t) 


(14) 


D(v,t) = 


00 

,,2 2 r 
16tt e 


m 


dk ««o ek - kv) £ o k (t) 


(15) 


T 2k (t) = f Ak 


3^ V v " t) 


CO 


V = 


ek 


(16) 


where <^(t) is the spectral energy density |E. | 2 /8 tt at the wavenumber 
k . A detailed discussion of the above system of equations can be found in 
Davidson (1972) . We give here only the features relevant to the problem 
under consideration. The unstable waves will have frequency and will 

be generated only in regions where Y2j c ( t ) > 0 > which corresponds to phase 


velocities in the region (Fig. 1) 


V, - AV, < V < V. . (17) 

b b p b 

At each instant of time half of the directed energy of the beam 
goes into electrostatic waves, and half into flattening the beam, i.e., 

<§(t) = jdk ^(t) * n b mV b AV fc (t) . (18) 

The beam will flatten as in Fig. 2 within a time of the order OWn^)^ 1 , 
after which all wave generation by the beam will cease, and the beam 
spread AV^ - . Applying quasilinear theory to the case 

of the Type III bursts we find that the electron beam decelerates without 
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a distance of ^30 km. Kaplan and Tsytovich (1968) and Smith and Fung 
(1971) discussed the possibility that induced nonlinear scattering of plasma 
waves out of the beam direction would stabilize the beam. This process 
was shown to be inadequate for electron beams (Zheleznyakov and Zaitsev 
1970a, Smith and Fung 1971). 

Let us examine now the nonlinear effects introduced by the 
oscillating two-stream instability. As the unstable waves grow in the 
region resonant with the beam, the background plasma sees at time t a 

ik 

field oscillating with frequency ^ and amplitude E Q (t) . When this 
field exceeds the threshold field given by Eq. (8), the oscillating two- 
stream instability starts pumping wave energy out of the mode k , into 
plasma oscillations with lower phase velocity (Fig. 3) and zero frequency 
(i.e., purely growing) ion density fluctuations. This effect removes plasma 
waves from resonance with the beam, and may cut off formation of the plateau. 
This will occur if the growth time of the new waves becomes shorter than the 
generation time of the resonant waves. 

We can now calculate the conditions necessary for stabilizing the 
beam against quasilinear diffusion. Our calculation is similar to that of 
Tsytovitch and Shapiro (1965) , and the reader is referred to that reference 
for more details. The procedure is straightforward. One finds an approxi- 
mate solution to the coupled equations for the beam-plasma and oscillating 
two-stream instabilities. This solution permits an estimate of the variation 

of the beam distribution function <5F. . 

b 

The beam distribution function evolves according to Eqs. (13) 
and (15). In the resonant region, the evolution of the field energy is 
described by 




bk 


2Y 2k ^ok 




9t 


(19) 
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where ; *- s t ^ e energy excited in the nonresonant region by 

the oscillating two-stream instability. In the nonresonant region the 
field energy is given by 

at 2y ik ^Nk 


( 20 ) 


Although we are interested in the early stages in the develop- 
ment of the instabilities, we will use Eq. (11) for y • This is because 

the growth rate at marginal stability is smaller than that given by Eq. 

2 2 

(10) well above threshold. In Eq. (11) we neglect E^ compared with E^ . 
[The validity of this assumption can be confirmed from the computer simu- 
lations of Kainer, Dawson and Coffey (1972)]. The resulting set of equa- 
tions (13) , (15) , (19) and (20) have the identical form to those solved 
by Tsytovich and Shapiro (1965). Consequently, the variation of F^ is 
given by 


3(6F ) 

p 

D 

9v 

- P&n 


^ ^ m V b A \ 

^(t = °) 


fb 

9v 


( 21 ) 


where fc^(t=0) is the initial energy density of the waves in the non- 
resonant region and p = y(k) = Y 2 k^ Y lk * 

The beam will stabilize against plateau formation if 


96F t 

9v 


« 


9v 


which implies that 


... r-b > i 

J 


« i 


( 22 ) 


Because the logarithmic dependence is very weak, we can replace it with 


/ 
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an upper bound, namely. 


£n 


yn_ m V, AV, 
b b b 

<%(t = 0) 


. f n b m V b ,3 
« tol-. ^ t - n e X D 

e e 


« in 


2 3 

VT n V 
b e b 


- 30 , 




where = KT /m and we have used (V, /V ) - 20, n - 10^ cm ^ 

e e be e 


a) - -10 
e 


8 


Therefore, we may approximate (22) by 


and 


y < 0.03 


(23) 


Because y depends only on n^/n^ anc * A V^/V^ we can find the value of 
AV^/V^ necessary for stabilization as a function of n j-/ n e • From Eq. 
(16) we find that . 


'2k 



(M 2 

U ) 

\AV./ ’ 

e 

b 

(ID 

and replace 

If we 

also assume 


(24) 


2 

In order to find y^ , we use Eq. (11) and replace E o /8tt by the value 

2 

n^ m V^AV^ as given by Eq. (18) 


then 


and 



(25) 


(26) 


This constitutes an extremely conservative estimate, because the instability 
threshold can be easily exceeded and the more appropriate growth rate will be 
that given by Eq. (10) . In that case the beam is stable as long as 
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^ > 33 ( — 

Vk W “ I 


2/7 


(27) 


8 10 

In Fig. 4, using - 5 x 10 cm/s, - 10 cm/s, we have plotted 

AV^/V^ as a function of » assuming that equality holds in Eq. (27). 

From Fig. 4 we can see that only if the ratio of the beam-to- 

-4 

plasma density is much larger than 10 is the oscillating two-stream 

instability not capable of stabilizing the beam against formation of a 
quasilinear plateau. However, for smaller density ratios the quasilinear 
velocity diffusion of the beam will cease early in the process and plasma 
wave energy will be concentrated in the regions of phase space shown in 
Fig. 3. 


IV. Discussion 

A detailed calculation of the radiation spectra produced by the 
oscillating two-stream instability will be presented in a future publication. 
Here we simply demonstrate that the observed power can be accounted for by 
this theory. For brevity we concentrate on the 2u> e radiation but similar 
considerations also apply to radiation at the fundamental. 

As described previously, the oscillating two stream instability 
will create symmetric [i.e., 4^(k, w ) «^(-k, ] oscillation spectra 

in the non-resonant region (Fig. 3), predominantly in the direction of the 
beam. These electrostatic oscillations can collide head-on with each 
other and radiate a transverse wave at 2iO e (Fig. 5) . The process is 
greatly enhanced by the fact that the electrostatic spectrum is an even 
function of the phase velocity. In contrast, the linear beam-plasma in- 
stability creates an asymmetric longitudinal spectrum, with phase 
velocities in the direction of the beam velocity. Subsequent emission at the 
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second harmonic requires intermediate scattering to produce plasma waves 
with phase velocities antiparallel to the beam (Smith 1970) . 

The energy density of plasma waves, W , produced by the oscil- 
lating two-stream instability will be determined by the mechanism that 
stabilizes it. In the discussion below we consider two possibilities. 

The first is that the oscillating two-stream instability stabilizes at 
field energy densities close to the threshold given by eq. (8). At thresh- 
old, y given by Eq. (11) is near zero and the system is stable. The 
second possibility for stabilization arises because early in the develop- 
ment of the oscillating two-stream instability the electric fields are 
well above threshold. These plasma waves will be Landau-damped by particles 
in the thermal plasma. This results in the formation of suprathermal tails, 
which have been observed in the computer experiments of Kainer, Dawson, 
and Coffey (1972) . The instability will stabilize when the Landau damping 
due to the tails becomes equal to or larger than the growth due to the oscil- 
lating two-stream instability. A third possibility, which we shall not 
consider further, is that resonance broadening effects characteristic of 
strong turbulence can lead to stabilization. An analysis of the discussion 
by Bezzerides and Weinstock (1972) indicates that one must extract a large 
portion of the beam energy in order for resonance broadening to be important. 
We do not find it necessary to do this and, consequently, shall not consider 
effects of strong turbulence any further. 

In general, regardless of the details of the stabilization mechanism 
for the oscillating two-stream instability, the electromagnetic radiation 
emitted at 20)^ is given by (Sturrock, Ball and Baldwin 1965) 


J(^,,U)) 'u 2u) ) 
- T e 


✓37 r 2 e 2 

2 3 
m c 


dk 






l- 


.2 , 2 

k T k 




-k)C^(k) 


(28) 
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where the radiated wave number is given by k^ ^ /3 (o^/c and is much 
smaller than the excited electrostatic wave numbers. Therefore Eq. 
(28) can be written as 


JC^j, , (0 


'v* 2(0 ) = 

e 


4/3 It e 


2 3 k 2 cos 2 0 sin 2 6 dk <^(k) ^(-k) . 


m c 


( 29 ) 


One can estimate the emission by taking the spectrum to be almost one 
dimensional with a peak at the wavenumber with the maximum growth rate, 
1^ = akp , where a << 1 and k^ = X ^ . Using ^(-k) - <^(k) and 
noting that 


W = dk ^ N (k) 


(30) 


Eq. (29) gives 


J(k , (0 ^ 2(0 ) 4/3 TT e k T w 2 ? 0 2 n 

T e = 2 — 3 — — 3 sin 9 cos 6 

m c Otk^ 

2 


2 3Ai(M nKT u 

a ' c / \ n e KT g / e e i 


sin 2 0 cos 2 0 . 


(31) 


The flux at the Earth is 


2ir V . r J(k_ , 2(0 )"1 

rad v ~ T e ' erg 

2 2(0 2 
R L e -* cinsec Hz ster 


(32) 


where V ^ is the radiating volume and R the distance of the source 

from the Earth. We have assumed that there is no damping of the radiation. 

7 “3 8 -i 

Using as typical parameters n g K 8 x 10 cm , % 5 x 10 s , 

V 'v 5 x 10^ cm s ^ , T 'o 2 x 10^ k > V , 'u 10^® crn"^ , R 'v 2 x 10^ cm, 
e e ^ rad * * 

a .1 , and defining 



- 18 - 


n KT 
e e 

-9 2 W 

we find F ^ 2 x 10 A — r — at the 80 MHz level. 

m Hz 

We should emphasize that we use these parameters only for illustrative 
purposes. A later, more detailed analysis will include a range of para- 
meters extending from less than a solar radius to 1 A.U. 

-21 -17 2 

Observations give F ^ 10 -10 W/m Hz . Therefore in order 

-6 -4 

to account for the observed flux we need X ^ 10 - 10 . Because 

-3 -2 

n^m^AV^/n^KT^ ^ 10 - 10 , only a small energy transfer is required 

to account for the observed power level. Therefore, it does not appear 
difficult to stabilize electron beams against quasi-linear diffusion, and 
at the same time to produce sufficient radiation to account for the observed 
Type III burst intensities. 

The determination of the value of A and the details of the radiation 
spectrum depend on the detailed nonlinear stabilization mechanism of the os- 
cillating two-stream instability, which will be discussed in a future publication. 
However some further estimates are possible. If electrostatic fluctuations 
with A 'V 10 8 -10 ^ exist in the non resonance part of the spectrum, higher- 
energy electrons in the ambient thermal distribution will interact with them 
via Landau damping and form long suprathermal tails. The time scale for this 
interaction is determined by the diffusion coefficient and is approximately 
given at the 80 MHz level by 

V 2 A V 

T ~ ^ P h ^ 5 x 10“ 8 A -1 s 

orD oru) Av , 

e ph 



where AV , ^ V . is the range of phase velocities over which the oscil- 
pn ph 

lating two-stream instability is unstable. Therefore for values of 
-6-4 

A 'Vi 10 -10 the plasma can form non-Maxwellian tails within times 

shorter than 1 millisecond. One then finds the marginally 

stable situation shown in Fig. 6. Marginal stability analysis (Papadopoulos 
and Pongratz 1973) shows that the tails have a velocity dependence 'v v n , 
with n ^ 1-3 . The existence of the tails will stabilize the oscillating 
two-stream instability, while the enhanced electrostatic fields in the region of 
the tails will prevent the re-development of the beam-plasma instability. In this waj 
the role of the beam will be to form and sustain these power law tails in 
the coronal plasma, accompanied by an enhanced level of electrostatic fluc- 
tuations. 

The level of fluctuations will be determined by the maximum of 
either the non-thermal equilibrium level due to the long-tails, or the 
level of fields in the non-resonant region necessary to prevent the beam 
plasma interaction from resuming. The value of A due to non-thermal tails is 



where A fc is the level of thermal noise. The approximate equality results 
from the formation of a high temperature Maxwellian tail. However, a power 
law such as results from marginal stability analysis will produce an even 
greater enhancement. The value A^ to prevent the beam-plasma instability 
from resuming is given by 
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where y = Y^/Y^ from Eqs. (24) and (11) and E T is given by Eq. (8) 
with now determined by the tails. The value of X = max (X^ , X ) . 

Notice that the value of X q ^ 10 ^ which is of the order required to 
explain the observations. 

Similar considerations can be applied with respect to radiation 

at 0) . 

e 

V. Summary 

We have demonstrated that the oscillating two-stream instability, 
induced by the fields created from the interaction of the electron stream 
with the coronal plasma, can stabilize the beam against quasilinear dif- 
fusion. This interaction creates a marginally stable situation in which the 
beam propagates essentially unaffected while sustaining in the ambient 
plasma suprathermal tails and suprathermal enhanced plasma oscillations. 
These oscillations can then scatter from ambient density fluctuations, 
and from each other, producing the observed radio emission at the plasma 
frequency and its second harmonic. 
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FOOTNOTE 

One may easily verify that for the parameters under consideration, 
the condition for a monochromatic pump spectrum 

<WYi - 6K(t)/Y 1 < 1 

is easily satisfied within a couple of growth times of the beam-plasma 
instability. 
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FIGURE CAPTIONS 

1. Region of unstable wave generation in the beam plasma interaction. 

2. Final form of the electron distribution function predicted by quasi- 
linear theory. Unstable wave generation has ceased. 

3. Band of phase velocities excited by the oscillating two-stream in- 
stability. 

4. Stability criterion of an electron beam as a function of AV, /V, and 

b b 

8 10 

n^/n e> using V g = 5 x 10 cm/s and = 10 cm/s. 

5. Scattering of two electrostatic plasma waves at frequency co^ radiating 
an electromagnetic wave at frequency 2a> e . 

6. Final form of the electron distribution function predicted by the para- 
metric oscillating two-stream instability. Enhanced emission continues. 
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